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Abstract 

We investigate the ground-state property of an e g -orbital Hubbard model at quarter filling on a zigzag chain by 
exploiting the density matrix renormalization group method. When two orbitals are degenerate, the zigzag chain 
is decoupled to a doble-chain spin system to suppress the spin frustration due to the spatial anisotropy of the 
occupied orbital. On the other hand, when the level splitting is increased and the orbital anisotropy disappears, a 
characteristic change in the spin incommnsurability is observed due to the revival of the spin frustration. 
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The magnetic property of frustrated spin systems 
has been one of the central issues for many years in the 
research field of condensed matter physics [1]. For ex- 
ample, it is well known that in the antiferromagnetic 
(AFM) Ising model on a triangular lattice, there oc- 
curs macroscopic degeneracy for possible spin config- 
urations in the ground state [2]. In general, however, 
such high degeneracy is lifted to suppress the spin frus- 
tration, since the lattice is deformed to lower the lattice 
symmetry due to the spin-lattice coupling. 

On the other hand, recently there has been a rapid 
increase of interest in the effect of the interplay of spin 
and orbital degrees of freedom [3]. It has been empha- 
sized that the orbital anisotropy plays a significant role 
to cause a variety of cooperative phenomena in realis- 
tic materials. In particular, in geometrically frustrated 
lattices, it is expected that orbital ordering occurs to 
affect the spin frustration due to the spin-orbital cou- 
pling, since the orbital anisotropy leads to the non- 
uniform exchange interactions. 

In this paper, to clarify the key role of the orbital 
anisotropy in geometrically frustrated lattices, we in- 
vestigate an e g -orbital Hubbard model on a zigzag 
chain with one electron per site (quarter filling). When 
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the Hund's rule coupling J is small, the ground state is 
found to be paramagnetic (PM) [4,5], which is relevant 
to a geometrically frustrated antiferromagnet. Here we 
study the property of the PM phase and set J=0 for 
simplicity. The effect of J has been investigated for an 
e g -orbital degenerate model [5] . 

The Hamiltonian considered here is given by 

H = Yl i 77'4 7 <r d i + a 7 ' CT - (A/2) J^ipia ~ pib) 
1,0.,^,^' ,cr i 

+ U Y1 PhtPlfl + U ' PiaPib, (1) 

i,7 i 

where dw(A&cr) is the annihilation operator for an 
electron with spin a in the 3z 2 — r 2 (x 2 — y 2 ) orbital at 
site i, Pi 7 <T=d? 7(T di 7CT , and Pi 7 =I] CT pi ra .t*^, is nearest- 
neighbor hopping between 7 and 7' orbitals along the a 
direction. Note that the zigzag chain is considered as a 
double chain connected by a zigzag path. The hopping 
amplitudes are given by t* a =t/4, t^ b =tf a = — v3t/4, 
t* b =3t/A for the double-chain direction and f" a =i/4, 
tab=t ba =y/3t/8, t bb =3t/16 along the zigzag path. Here- 
after, t is taken as the energy unit. A is the level split- 
ting between 3z 2 — r 2 and x 2 —y 2 orbitals, U is the in- 
traorbital Coulomb interaction, and U' is the interor- 
bital Coulomb interaction. 
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Fig. 1. Orbital structure for (a) A— and (b) A — 1 from the 
DMRG results for JV=40 at U' = U=20. Spin correlation for 
(c) A— and (d) A— 1. Electron densities in (c) 3z 2 — r 2 and 
x 2 — r 2 orbitals and (f) optimal a and b orbitals. 

We analyze the model with N sites in the open 
boundary condition by using the density matrix renor- 
malization group (DMRG) method [6]. To reduce the 
size of the superblock Hilbert space, we treat each or- 
bital as a site. We employ the finite-system algorithm 
with keeping up to 200 states per block and the trun- 
cation error is estimated to be 4xl0~ 6 at most. 

In order to determine the orbital structure, we in- 
troduce new operators by using an angle 9i such as 
di d a = e i9i/2 [cos(f9 i /2)d iaCT +sin(f5 i /2)d if)CT ] and d iia = 
e Wi/2 [- sin(f3i/2)d iaCT + cos(f9i/2)dii, CT ] [7]. The optimal 
orbitals, a and b, are determined so as to maximize 
the orbital correlation T(q) = ]T\ j e ,q ' (i - j) {f?f?)/N 2 

withfT = £ CT (4.^-^4 j/2- In Figs. 1(a) and 
(b), the orbital structure is shown for A=0 and 1, re- 
spectively. When two orbitals are degenerate for A=0, 
orbital degree of freedom is active, but a 3a; 2 — r 2 orbital 
is selectively occupied to suppress the spin frustration, 
as shown in Fig. 1(a). Namely, the orbital shape ex- 
tends just along the double-chain direction, and the 
zigzag chain is decoupled to a double-chain spin sys- 
tem due to the orbital anisotropy. In fact, the ratio of 
the AFM exchange interaction along the double-chain 
direction J2 to that along the zigzag path Ji is esti- 
mated as J2 / Ji=64: 2 . On the other hand, for A=l, a 
lower-energy 3z —r orbital is favorably occupied, as 
shown in Fig. 1(b). Note that when the 3z 2 — r 2 orbital 
is fully occupied for infinite A, the orbital anisotropy 
disappears in the xy plane, i.e., J2/J1 — I, and the spin 
frustration becomes effective. 

In accordance with the variation in the orbital shape, 



the spin state is also changed. To clarify this point, it is 
convenient to reduce the present model to a spin system 
on the orbital-ordered background. Then, the present 
system is described by the zigzag spin chain, in which 
the spin correlation has a commensurate peak at q=n 
for 0< J2/ Ji<l/2, but the peak is gradually changed to 
an incommensurate one for J2/Ji>l/2, and eventually, 
we find the incommensurate peak at q—n/2 for infinite 
J2/J1 [8,9]. In Fig. 1(c), we show our DMRG result 
of the spin correlation S(q) = £ tJ e <q - (i - j) (5f Sf)/iV 2 
with Si = J2-,(pi~il ~ Pi7l)/2 for A=0. We find a clear 
peak at q=n/2, consistent with that of the zigzag spin 
chain with large J2/J1, since J2/Ji=64 2 for A=0. On 
the other hand, as shown in Fig. 1(d) , the peak position 
for A=l changes from q=ir/2 toward q=n, expected by 
analogy with the zigzag spin chain, since we estimate 
J2/Ji = 1.61 for A=l. The detail of the A dependence 
will be discussed elsewhere in future. 

Finally, let us consider how the orbital state changes 
in the intermediate region. In Fig. 1(e), we show the A 
dependence of the electron densities n 7 = J2i(pi~/} /N 
in 3z 2 — r 2 and x 2 — r 2 orbitals. With increasing A, elec- 
trons are forced to accommodate in the lower 3z 2 — r 2 
level, but the electron density in each orbital is found 
to change gradually without any singularity. To under- 
stand this behavior, we evaluate the electron densities 
for the optimal orbitals. As shown in Fig. 1(f), it is 
found that one of the optimal orbitals is occupied irre- 
spective of A and the fluctuation is very small even in 
the intermediate region. Namely, the present system is 
always regarded as a one-orbital system, although we 
have considered the multi-orbital system. 

In summary, for A=0, the 3a; 2 — r 2 orbital is selec- 
tively occupied to suppress the spin frustration and the 
zigzag chain is decoupled to a double chain due to the 
orbital anisotropy. For large A, the 3z 2 — r 2 orbital is 
occupied and the spin frustration revives, leading to 
the change in the spin commensurability. 
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